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Abstract 

D 

Ph . 

Using a 73.6 pb ^ data sample of T(2S') events collected with the CLEO 
I II detector at the Cornell Electron Storage Ring, we have investigated the 

hadronic transitions between the T(25) and the T{1S). The dipion tran- 
■ sition T(25) — > T(lS')7r''"7r~ was studied using two different analysis tech- 

'sj" ■ niques. Selecting events in which T(15) — > e~ , iJ,~ ("exclusive" anal- 

, ysis), and using the T(IS') leptonic branching ratios world averages from 

(M ' the PDG review, we found ^(T(25) ^ T(15)7r+7r-) = 0.189 ± 0.004 ± 0.010, 

O 

00 



while using a method allowing T(15) — > anything ("inclusive" analysis) 
0> ■ we found ^(T(25) ^ T(15)7r+7r^) = 0.196 ± 0.002 ± 0.010. The appro- 

■ priate average of the two measurements gives B{T{2S) T(15)7r"'"7r~) = 

^ . 0.192 lb 0.002 lb 0.010. Combining the exclusive and inclusive results we de- 

Ph! rive the T(15') leptonic branching ratios Bee = 0.0229 ib 0.0008 ib 0.0011 and 

^ ; Sf^f^ = 0.0249 lb 0.0008 ib 0.0013. We also studied T(25) ^ T(15)7r07r0 and 

obtained i3(T(25) T(lS')7r°7r°) = 0.092 ib 0.006 ib 0.008. Parameters of the 
■H ■ TTTT system (dipion invariant mass spectra, angular distributions) were ana- 

. lyzed and found to be consistent with current theoretical models. Lastly, we 

5^ I searched for the ry and single vr*^ transitions and obtained the upper limits 

B{T{2S) T{lS)ri) < 0.0028 and B{T{2S) T{1S)tt°) < 0.0011. 
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I. INTRODUCTION 



The hadronic transitions in heavy quarkonium systems provide an experimental testing 
ground for the theoretical calculations of nonperturbative QCD [|I| and can give informa- 
tion on the structure of QCD confinement as well as on the gluon content of light hadrons. 
Historically, studies of the hadronic transitions T(25') T(lS')7r7r were preceded by investi- 
gations of the transitions r( — > 777r7r and ^z^' — ^ ^/'tttt. All three are examples of A/ = dipion 
transitions. In the decay Vj — > t^tttt the pions fit reasonably well to a phase space mass 
spectrum 0. Soon after the discovery of charmonium Q], and the subsequent observation 
of the il)' ip-KiT transition, it was found that in this transition the dipion invariant mass 
spectrum cannot be adequately described by a phase space mass spectrum. The challenge of 
providing an acceptable description of the observed data attracted considerable theoretical 
attention. With the discovery of another family of heavy quarkonium states, the family of 
T resonances, the theoretical calculations were extended to include bottomonium. 

Figure |1] shows the bottomonium levels up to the T(25') and possible transitions between 
them, including radiative and rare (Svr and single ir^) transitions |Q. The hadronic transitions 
between the bottomonium levels are soft processes (typical transition energies are 0.3-0.9 
GeV) and are thereby difficult to treat perturbatively. Typically, the heavy quarkonium 
hadronic transition (qq)' — {qq)X is treated as the factorizable product of two processes: 
first, the transition from (qq)' to (gg) with the emission of gluons (usually two), followed by 
the hadronization of the gluons to the state X (i.e., the production of X from the vacuum 
in the presence of the gluon color field). 




FIG. 1. Transitions in the bottomonium. 

Although nonperturbative, the hadronic transitions between heavy quarkonia can nev- 
ertheless be described in the context of a "multipole" expansion scheme where the gluon 
fields are expanded in a multipole series, similar to the electromagnetic transitions, as first 
outlined by Gottfried [^. In the framework of the multipole expansion, Yan 0, and later 
Zhou and Kuang ^ calculated the transition rates and derived a parameterization for the 
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dipion invariant mass spectrum in the T{2S) —>■ T{1S)tt7[ transitions. They used the quark- 
confining string model |@] to describe the intermediate state of the hadronic transition and 
calculate the hadronization matrix element. Rather then writing the gluonic degrees of free- 
dom for the quark-confining string, Voloshin and Zakharov (VZ), and afterwards in a 
revised analysis Novikov and Shifman (NS), used an alternate approach and wrote the 
general form of the QCD field tensor in the chiral limit to obtain the hadronization matrix 
element. In both approaches the hadronization matrix element is constrained by current al- 
gebra, partial conservation of the axial current (PCAC), and gauge invariance. The essential 
mass dependence of the matrix element is very similar in all cases: it vanishes for dipion 
mass approaching threshold, and peaks at larger values of m^rTr. In the NS and VZ models, 
the model parameters are derived from "first-principles" , as opposed to the Yan et al. model 
where the parameters are determined phenomenologically from a fit to ip' — > ipiiTc. 

The results presented in this paper were obtained using the world's largest available 
data sample of T{2S) decays (73.6 pb~^ of integrated luminosity on-resonance, and 5.2 pb~^ 
off-resonance) collected with the CLEO II detector at the Cornell Electron Storage Ring 
operating at the T{2S) center of mass energy in December 1994. Similar investigations 
were performed by several collaborations including ARGUS [0, CUSB |T^, CLEO |jl3[ 



and Crystal Ball |T^. Our data sample is larger by at least a factor of two in integrated 
luminosity than each of the previous measurements, with the number of T{2S) resonant 



decays Nr(2S) = (488 ± 18) ■ 10^ p 



II. DETECTOR 



CLEO II is a general purpose detector 16 for measuring charged and neutral particles 



in the energy range from ^ 50 MeV to ~ 6 GeV. Its three concentric wire drift chambers, 
covering 95% of the sohd angle, detect charged particles and perform particle identifica- 
tion using specific ionization energy loss measurements (dE/dx) in the outer chamber. A 
superconducting coil provides a magnetic field of 1.5 Tesla, giving a momentum resolution 
of {Sp/pY — (0.0015p)^ + (0.005)^, where p is the momentum in GeV/c. A time-of-fiight 
system, just outside the drift chambers, consists of plastic scintillation counters and serves 
as a primary triggering system; it also provides some particle identification information. Be- 
yond the time-of-fiight system, but inside the solenoid, is an electromagnetic calorimeter, 
consisting of 7800 thallium-doped Csl crystals arranged as two endcaps and a barrel region. 
The central barrel region of the calorimeter covers 75% of the solid angle and achieves an 
energy resolution of 6E/E{%) = 0.35/i?*^'''^ + 1.9 — O.IE, where E is the shower energy in 
GeV. The endcaps of the calorimeter extend the sohd angle coverage to about 95% of An, 
although energy resolution is not quite as good as in the barrel. Proportional tracking cham- 
bers for muon detection are located in between and outside of the iron slabs that provide 
the magnetic field fiux return. 

In our analysis we used a customized version of LUND/JETSET [|l^] program as a Monte 
Carlo event generator. The final state particles are then propagated through and decayed in 
the CLEO II detector using a GEANT [jl8i based detector simulation. 
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III. TRANSITION T(2S') ^ T(15)7r+7r- 



We studied the dipion transitions T{2S) T{lS)TnT using two different techniques. The 
first one selects events where there is a e^e^ or fi^fj^^ pair recoihng against the dipion system, 
which is assumed to result from T(15') e^e~,fi~^fi~ ("exclusive" measurement). In the 
second technique we do not specify how the T(IS') decays, selecting all events which have 
a TT+TT^ pair ("inclusive" measurement). The two measurements are complementary to each 
other and provide us with important cross-checks. 

A. Exclusive final states with T(1S) e^e^ 

We use the following selection criteria for the exclusive events with tt+tt^/^/^ in the 
final state. We demand four tracks in the event which pass track quality requirements, two 
of them (the lepton candidate tracks) must have momentum greater than 3.5 GeV/c and 
come from the vicinity of the interaction point which is defined as a cylinder of 3 mm x 
10 cm (radius x length) with its axis lying along the beam direction, and the other two 
(pion candidate tracks) must have momentum less than 0.5 GeV/c and come from a similar 
cylindrical region 4 mm x 12 cm around the interaction point. To suppress background from 
radiative Bhabha events with 7-conversion we require that the cosine of the angle between 
the pion tracks satisfy cos6'7r7r < 0.9. We identify electrons by the combined requirement 
that the ratio of the shower energy to the momentum of the matching track is close to 
1 and that the lateral energy deposition in the calorimeter is consistent with the electron 
hypothesis. Events with muons are identified by requiring that the sum of the maximum 
penetration depths of the two tracks into the muon system absorber be greater than four 
hadronic absorption lengths.[| 

The missing mass Mjniss = \J {Mr (2s) — -^^tttt)^ — pI^ (i-^-, the mass recoiling against the 
dipion system) distributions for both ee and /i/i channels are shown in Figure 0. One can 
see a clean signal with very little background in the side-bands,0 thus we use a simple event 
count to obtain the number of observed events both in Monte Carlo (to calculate efficiencies) 
and in data. 

The three largest sources of background are QED radiative processes with 7-conversion, 
two-photon double-tag production of vrvr (in the ee channel) and one-prong r-decays from 
T(IS') TT. Due to our high lepton momentum and lepton identification requirements the 
contamination from r-decays to our data sample (which we directly subtract from the number 
of observed events) is very small: less then one event in each channel considered. To eliminate 
QED radiative and two-photon background we use the method of side-band subtraction: we 
count the number of events in the side-bands of our signal region and extrapolate this number 



^As one can notice we just require for two high momentum lepton tracks to be present in the event 
and do not apply any further criteria on the dilepton system. 

^The signal region is defined as the missing mass interval (9.43,9.49) GeV, the side-bands are 
defined as (9.20, 9.40) U (9.52, 9.70) GeV in both dilepton channels. 
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FIG. 2. The missing mass distributions in the exclusive T(2S) T(lS')7r+7r measurement. 

into the signal region. In this way, we find the background contamination to be 8.7 events 
(0.9%) in the ee channel and 3.8 events (0.3%) in the /i/i channel. 



TABLE I. Numbers of events observed after background subtraction, efficiencies, prod- 
uct of branching ratios B(T{2S) T(15)7r"'"7r~) • B{T{1S) and branching ratio 
i3(T(2S') — > T(15)7r+7r~) for the exclusive measurement. 





j\j-observed 




B^^-Bu (xlO-3) 


B 


ee 


956.2 ± 30.9 


43.7 ± 1.4 


4.5 ±0.1 ±0.2 


0.178 ±0.006 ±0.015 




1130.1 ±33.6 


47.5 ± 1.6 


4.9 ±0.1 ±0.2 


0.196 ±0.006 ±0.011 



Knowing the efficiencies eu from the Monte Carlo simulation,^ we can calculate 
the products of two branching ratios S(T(25) ^ T(lS)7r+7r-) ■ B{T{1S) l+l') = 
A^,f^"™^'^/(enA^T(25)), as shown in Table |. Using the Particle Data Group (PDG) values 



for B{T{1S) e+e") = 0.0252±0.0017 and B{T{1S) /i+zz") = 0.0248±0.0007 we obtain 
the T{2S) dipion branching fraction. Combining the results from both channels we find: 

B{T{2S) T(15)7r+7r-) = 0.189 ± 0.004 ± 0.010 

where the first error is statistical and the second is systematic]^ (for a breakdown of systematic 
errors see Sec.V). 



^For all our sub-analyses we used the Voloshin and Zakharov ^ model with A = 3.44 to generate 
the dipion invariant mass spectrum in the Monte Carlo simulation. 

^When we average over the two dilepton channels, we treat correlated and uncorrelated errors 
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In Table ^ we compare our result with other exclusive measurements. 

TABLE II. B{T{2S) T(lS)7r+7r-) • e(T(lS) ^ l+r ) in units of IQ-^. 



ARGUS n 



Crystal Ball 14] 



CUSB H 
CLEO ll3| 
LENA IT 



this analysis 



4.4 ±0.2 ±0.4 
4.9 ±0.4 ±1.0 
5.4 ±0.3 ±0.4 
5.4 ±0.4 
6.1 ±2.3 

4.66 ±0.10 ±0.23 



average 



4.82 ±0.18 



B. Inclusive final states 

In our inclusive analysis of T(2S') —>■ T(lS')7r+7r^,T(15') ^ X we select events with at 
least two tracks which pass our track quality requirements, have momentum less than 0.5 
GeV/c and come from the interaction region. We also require that the invariant mass of the 
two pion candidates lie between 0.27 GeV/c^ and 0.57 GeV/c^. 




FIG. 3. Missing mass distribution from the inclusive T(25') T(15)7r"'"7r events: a) the full 
distribution, b) the region near the T(IS') mass, with the fit to the T(IS') peak. 

The signal appears in the missing mass plot shown in Figure ^ along with the fit to the 
T(IS') peak. The fitting function we use is a double-gaussian with the difference between 
the means fixed at zero, plus a third order polynomial for the background. The number of 



separately in calculating the overall systematic error. 
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fitted events in the peak is Ni^d = 50566 ± 570. The efficiency for the T{2S) T(lS')7r+7r^ 
channel has been calculated from a Monte Carlo simulation and determined to be tind = 
(52.9 ± 2.0)%. From these two numbers and the total number of T{2S) produced we find 
the branching ratio for the transition T(2S') — > T(lS')7r+7r~: 

B{T{2S) T(15)7r+7r-) = — — = 0.196 ± 0.002 ± 0.010 

A comparison of this result with previous inclusive measurements is given in Table 



TABLE III. B{T{2S) T(lS')7r+7r ) inclusive measurements. 



LENA PI 
ARGUS (ll 
CLEO in 
this analysis 



0.26 ±0.13 
0.181 ±0.005 ±0.010 
0.191 ±0.012 ±0.006 
0.196 ±0.002 ±0.010 



average 



0.190 ± 0.007 



Combining the results of the exclusive and inclusive measurements, and taking into ac- 
count correlations between the systematic errors, we obtain:^ 

B{T{2S) T(15)7r+7r-) = 0.192 ± 0.002 ± 0.010 

Knowing the number of inclusive and exclusive events, we can now calculate the T(15') 
leptonic branching ratios B{T{1S) l^l') = {Nueind) / {Nincieu): 

Bee = B{T{1S) e+e") = 0.0229 ± 0.0008 ± 0.0011 

B^^ = B{T{1S) = 0.0249 ± 0.0008 ± 0.0013 

which agree well with the corresponding PDG values. 

IV. TRANSITION T(25) T(15)7r°7r° 

To analyze the transition T(2S') — >■ T(15')7r°7r° exclusively in the final states with 
T(IS') — >• e"'"e~, /i"*"/!", we apply criteria similar to those used in our T(2S') — >• T(lS')7r+7r^ 
exclusive analysis: we require two good quality charged tracks passing the criteria for a 
lepton candidate, and the same e//i-identification criteria. 

We reconstruct tt" candidates from photon showers in the calorimeter. The photons are 
required to satisfy the following criteria: (1) the absolute value of the cosine of the polar 
angle (the angle between the photon and the beam directions) should be less than 0.95 to 



^This number is subject to the use of the PDG values for the T(15) leptonic branching ratios. 
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exclude the region of "hot" (noisy) crystals in the endcaps close to the beampipe, (2) the 
photon energy must he in the interval 0.05 GeV < < 0.43 GeV, (3) the angle to the 
closest projected charged track should be greater than 15°, (4) the shower should not be a 
fragment of a larger shower, and (5) the pattern of energy deposition should be consistent 
with the single photon hypothesis. Photons satisfying these requirements are combined into 
pairs to form vr'' candidates. Combinations with momentum greater than 0.5 GeV/c are 
excluded from further consideration. The pair of vr^'s remaining with the minimal value of 
the pull ^>S'^i72 '^7374' where S^^ = (m^^ — m^o)/(Tm^^ is then selected, and the missing 
mass calculated (Figure ||). As is the case with charged pions we see clean signals in both 
lepton channels. Because of the poorer momentum resolution of reconstructed tt'^'s than that 
of charged vr's, the distributions are considerably wider. 
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FIG. 4. The missing mass distributions in the exclusive T{2S) T(lS')7r vr measurement. 



TABLE IV. Numbers of events observed after background subtraction, efficiencies, prod- 
uct of branching ratios B{T{2S) — > T(lS')7r''7r'^) • ;B(T(1S') l^l^) and branching ratio 
B{T{2S) T(lS')7r''7r'^) for the exclusive measurement. 





j\j'observed 




B^^-Bu (xl0"3) 


B 




ee 


133.2 ± 11.5 


12.3 ± 1.0 


2.2 ±0.2 ±0.2 


0.088 ± 0.008 


±0.010 


nil 


142.5 ± 11.9 


12.2 ± 1.0 


2.4 ±0.2 ±0.2 


0.096 ± 0.008 


± 0.009 



Once again we perform a side-band subtractionFI to extract the number of observed events 



^Here the side-bands are (9.10,9.40) U (9.55,9.80) GeV in both channels and the signal region is 
(9.40,9.55) GeV. 
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(we estimate the background to be 3.8 events, or 2.9%, in the ee channel and 1.4 events, or 
1.0%, in the /i/z channel). 

The yields and efficiencies for exclusive T{2S) T(lS')7r'^7r°,T(lS') l^l^ transitions 
are presented in Table From these numbers we calculate the product of branching 

ratios i3(T (25) ^ T{1S)tt%°)-B{T{1S) /+/-) = A^^f ™7(eHArx(25)), and using the PDG 
values for T(IS') l^l^ , we find B{T{2S) T(lS')7r°7r°) which is also reported in Table |y. 
Averaging over the two dilepton channels, we obtain: 

B{T{2S) T(lS)7rV) = 0.092 ± 0.006 ± 0.008 

In Table previous determinations of B{T{2S) T(15)7r°7r°) ■ B{T{1S) /+/") are 
compared. From our two exclusive measurements we find the ratio B{T{2S) — >■ 
T(15)7r°7r0)/i3(T(25) ^ T{lS)TT+n-) = 0.49 ±0.06 which is close to the isospin zero expec- 
tation of 0.53. 



TABLE V. B{T{2S) T{1S)tt^tt^) ■ B{T{1S) l 


) in units of 10 ^. 


ARGUS 11 
Crystal Ball |14] 

cusB m 

this analysis 


2.3 ±0.4 ±0.5 
2.3 ±0.3 ±0.3 
2.9 ±0.5 ±0.3 
2.29 ±0.14 ±0.20 


average 


2.34 ±0.19 



An inclusive analysis of the T(2S') T(lS')7r°7r° transition gave a numerically consistent 
result, however because of the enormous combinatoric background, this measurement has 
very httle statistical weight. 



V. TRIGGER EFFICIENCY AND SYSTEMATIC ERRORS 

The trigger system of the CLEO II detector, described in detail elsewhere [pi 



designed for efficient triggering of two-photon, tau-pair, and hadronic events. There were 
eight active trigger lines in total during the T(2S') data taking, but only four of them 
are important in selecting events containing approximately back-to-back electron or muon 
pairs plus additional energy clusters in the calorimeter. To ffie, these trigger lines require 
either two hits in the opposite hemispheres in the time-of-fiight system or in the calorimeter, 
or a hit in the time-of-fiight barrel region plus a track in the vertex detector (with small 
variations from line to line). Our estimates of the overall trigger efficiencies from a Monte 



Carlo simulation of the trigger system are reported in Table ^ . 

TABLE VI. Trigger efficiencies. 





T(25) ^ T{IS)ti+ti- 




T{2S) 


T(15)7r07r0 




T(15)^ee T{IS) ^ 


T(15) ^ X 


T(15) ^ ee 


T{is) ^ 


Efficiency 


0.961 ± 0.008 0.962 ± 0.015 


0.990 ±0.011 


0.982 ±0.031 


0.977 ± 0.042 
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T'A RT T? \nj 
iAJjijll/ Vll. 


Sources and maj 


;nitudes of systematic 


errors. 






Systematic error (%) 




T{2S) 


T{IS)tt+tt- 




Source 


Exclusive 


Inclusive 


Exclusive 


Multiplicity of event 


— 


2.0 




Trigger efficiency 


0.9 / 1.6" 


1.1 


3.1 / 4.2 


Tracking 


2.8 


2.8 


— 


TT ''-finding 


— 




7.0 


Finite MC sample 


0.5 


0.5 


0.5 


Background subtraction 


0.3 / 0.2 




1.5 / 0.9 


Leptonic branching ratios 


6.7 / 2.8 




6.7 / 2.8 


Fitting function 




0.5 




ivrProd 


3.7 


3.7 


3.7 


Total 


8.2 / 5.7 


5.2 


10.9 / 9.4 



separately for ee/nfj, channels. 



The dominant systematic errors in our analysis come from uncertainties in the total num- 
ber of produced T{2S) resonance events, the leptonic branching fractions of the T(15'), and 
the charged track and 7r° finding efficiency. Other systematic errors are due to uncertainties 
in trigger efficiencies, event environment effects, background subtraction, and the shape of 
the fitting function (inclusive analysis only). The complete breakdown of systematic er- 
rors is given in Table |V11| (relative errors in percent). All these errors are considered to be 
uncorrelated and separately contribute to systematic uncertainties in our branching ratios. 



VI. DIPION INVARIANT MASS SPECTRA IN T(2S') ^ T(15)7r7r TRANSITIONS 



There have been several theoretical predictions for the dipion invariant mass distribution 
since a significant difference from phase space was found in ip' J/ipTrir transitions 



As shown in Figure ^, the dipion transition is treated as a factorizable two-step process: 
emission of gluons from heavy quarks and conversion of the gluons into light hadrons. The 
dipion invariant mass spectrum is determined by the second step, in the hadronization of 
two gluons emitted by the decaying bottomonium - a process which is not well understood. 



3031197-033 

Hadrons 




Emission 



FIG. 5. A hadronic transition as a two-step process. 
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The following parameterizations were used in fitting our experimental distributions: 
• Yan model: 



^^PS- [(mL - 2mlf + ^(mL - 2ml) " - Ami + 2K\l+^^)) + 0{^)] 



where K = ^^-S^- 
Voloshin and Zakharov model: 

da 



dm.,. 



oc PS ■ [ml 



Am 



2i2 



Novikov and Shifman ITQ] model: 
da 



dm,. 



2m^ 

oc PS ■ [ml^ - k{M2 - M,f{l - ^) + 0{k^)f 



In all the above formulas M2 = Mx(25), Mi = Mr (is) and PS is the phase space factor: 



PS 



[mt 



AmDlMf + Mi + mi^ - 2{Mfml^ + M^ml^ + MfM|; 



4M| 



A. The TT^TT invariant mass spectrum 

We extract a dipion invariant mass spectrum from both the inclusive and exclusive sam- 
ples. The dipion invariant mass spectrum from exclusive events is shown in Figure where 
we have combined results from both ee and nn channels. The data points in this histogram 
are the sideband-subtracted yields for the corresponding bins in m^^^^, where each data point 
has been corrected for acceptance (Figure H). The fits to this spectrum, using the afore- 
mentioned parameterizations are also shown in Figure they are all consistent with our 
data. 

To produce the dipion invariant mass spectrum in the inclusive measurement, we use 
a two-dimensional plot of m^^r vs. Mmiss (shown in Figure ^ which we slice in bins of 
mnn, project onto the Mmiss axis and then fit the projections with a double Gaussian for the 
T(15') peak plus a third order polynomial to represent the background. We correct the fitted 
number of T(IS') events for acceptance on a bin- by-bin basis (Figure ^) to obtain the dipion 
invariant mass spectrum. Fits to the the different parameterizations of this distribution are 
given in Figure In Table |VII1| we have compiled the values of the fitting parameters, their 
errors, and values of the fits for both the exclusive and inclusive measurements. 
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FIG. 7. Dipion invariant mass spectrum from exclusive T(25') — >■ T(lS')7r"'"7r events 
(corrected for acceptance). 

B. The TT^TT^ invariant mass spectrum 



We have a sample of T{2S) — ^ T(lS')7r°7r*^ exclusive events, large enough to produce the 
n^TT^ invariant mass spectrum for this transition. Here too the data points in our histogram 
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3031197-034 



- T(2S)-^T(1S)7r+7r-,T(1S)-^X 




FIG. 9. Dipion invariant mass spectrum from inclusive T{2S) — T(lS')7r"'"7r events 
(corrected for acceptance). 
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TABLE VIII. Fit results for the 7r"'"7r invariant mass spectrum parameterizations. 





Exclusive events 




Inclusive events 


Model 


Fit parameters 


/ 13 d.f. 


Fit parameters 


XV 13 d.f. 


Yan § 


B/A = -0.132 ±0.018 


15.6 


-0.154 ±0.014 


9.3 


Voloshin and Zakharov 


A = 3.11 ±0.18 


17.5 


3.42 ±0.16 


6.6 


Novikov and Shifman [10| 


k = 0.138 ±0.009 


15.1 


0.153 ±0.008 


8.4 



are the yields of exclusive /+/ tt^tt^ events in each corresponding m,r7r bin, corrected for 
acceptance (see Figure The fits to the acceptance-corrected tt^tt^ invariant mass spectrum 



TABLE IX. Fit results for the vr^vr'^ invariant mass spectrum parameterizations. 



Model 




Exclusive events 
Fit parameters 


XV 13 d.f. 


Yan § 

Voloshin and Zakharov 
Novikov and Shifman 


i 

10| 


B/A = -0.145 ±0.040 
A = 3.35 ±0.49 
A; = 0.139 ±0.022 


10.8 
11.1 
10.9 



are shown in Figure 0, with fit results reported in Table pC . 
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FIG. 10. Dipion invariant mass spectrum from exclusive T(2S') — > T(lS')7r''7r^ events 
(corrected for acceptance). 
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C. Combined results for the vrvr invariant mass measurements 



In order to compare the results of our analysis with the results of other experiments, we 
perform a simultaneous fit to the exclusive and inclusive 7i~^7t~ invariant mass spectra. We 
do not include the tt^h^ measurement in the combined fit because it has a slightly different 
parameterization (due to the difference in mass between neutral and charged pions) and much 
lower statistical significance. The fits to the combined data of the exclusive and inclusive 
T{2S) T(lS')7r+'7r~ decays are shown in Figure |Tl|. In Table |^ we compare the results of 
our combined fit with the results from previous experiments. 



TABLE X. Values of fit parameters using different parameterizations of the tttt invariant mass 
spectrum. 



Model 


Yan i 


Voloshin and Zakharov 


9 


Novikov and Shifman 


m 


Parameter 


B/A 


A 




k 





Crystal Ball 11 -0.18 ±0.15 3.3 ± 1.2 0.14 ±0.05 

CLEO m -0.18 ±0.06 3.2 ±0.4 0.15 ±0.02 

ARGUS -0.154 ±0.019 3.30 ±0.19 0.151 ±0.009 

this analysis -0.145 ±0.011 3.28 ±0.12 0.146 ±0.006 
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FIG. 11. Combined fit to the dipion invariant mass spectrum from exclusive and inclusive 
T{2S) T{IS)tt+tt- events. 
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VII. ANGULAR DISTRIBUTIONS 



The angular distributions in vrvr transitions were studied using our exclusive and inclusive 
TT+TT^ data samples. In e^e~ annihilation the T{2S) is produced polarized with its spin axis 
lying along the beam axis. This total angular momentum (and its projection onto the beam 
axis) must be conserved. There are three possible angular momenta in the final state of the 
dipion transition (Figure |12D: the total spin of the T(15') J, the internal orbital angular 
momentum of the dipion system 1 (the total spin of the dipion system s = 0) and the orbital 
angular momentum of the dipion system relative to the T(IS') L |24|. 




L = even 



Y(2S) 



FIG. 12. Angular momenta in the tttt transitions. 



Since the transition is expected to be dominated by El - El gluon radiation, the angular 
momentum of the hb system is not changed by the dipion decay and the polarization of the 
parent T{2S) should be observed in the subsequent decay of the daughter T(IS'). This is 
verified in the cos 9 and distributions of the outgoing with respect to the beam shown 
in Figure [l^: the expected (1 + cos^ 0) distribution is clearly seen in the cos^/+ distribution 
and the azimuthal distribution (f)i+ is reasonably fiat.[] 

The quantum numbers of both the T{2S) and T(15') are J^*" = 1 and J*^ = 0~; the 
dipion system has J*^*^ = 0^^. Parity forces / and L to be both even or both odd. The 
G-parity for the dipion system^ is 1 and from the formula G = (—1)'+''+^ with / = 0, s = 
0, G = 1 we find that /, hence L, must be even. 

All theoretical models describing the dipion invariant mass spectrum predict the pions to 
be emitted predominantly in the s-wave state (/ = 0), although there exists a prediction for 
the d— wave contribution (/ = 2) |jlO[ of the order of 1%. The c?— wave contribution can be 
observed in the cos 6* distribution of the angle 6* of the vr"*" in the tttt center of mass frame 
with respect to the tttt direction. (See Figure |T3| for definitions of angles.) This is shown 
in Figure |l^ along with the (p* distribution which should be fiat. It is possible to fit the 



^It should be flat because CESR beams do not have azimuthal polarization. 



®The operation of charge conjugation followed by isospin rotation does not change the state of the 
dipion system. 
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FIG. 13. Frames of reference and definitions of angles for the vrvr transitions. 



3031097-024 




FIG. 14. Angular distributions of from T(IS') l^l in the center of mass frame (corrected 
for acceptance). Solid lines are dN/d{cos 9i) = N ■ {1 + cos^ 9i) and dN/d(j)i = const fits. 



cos 6** distribution for our exclusive data sample to a coherent sum of s- 
here represents the size of the c?— wave contribution: 



and c?— waves; e 



dN 



d{cosei] 



oc 



0|2 



with the fit result: e = 0.077 ± 0.041. We studied whether this effect is present in our 
inclusive T{2S) T{1S)tt^tt^ sample (Figure |16|a). In the inclusive measurement the fit 
result is: e = 0.028 ± 0.027. Performing a simultaneous fit to the combined data from the 
exclusive and inclusive measurements (Figure ^6|3) we find: 



0.042 ± 0.022 
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0.15 



3031197-025 




COS 

FIG. 15. cos 9* and i^* distributions of tt"*" in the center of mass frame of tt~^tt~ system 
in the exclusive T{2S) — > T(15)7r"'"7r~ measurement (corrected for acceptance). Sohd hnes are 
dN/d{cos 9*) = N- \^A^^Y^ + eY^l"^ and dN/d^* = const fits. 
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FIG. 16. a) Fit to the cos6* distribution in the inclusive T(25) — > T(lS')7r^"7r measurement, 
b) Combined fit to the cos 6* distributions (corrected for acceptance) . 



Our results demonstrate the strong s— wave dominance expected in the dipion transition and 
show some indication of a (i— wave contribution on the order of a few percent. In a similar 
analysis, ARGUS ^ obtained e = 0.018l[i:J|]^. 

To examine further the question of a possible (i— wave contribution we performed a fit to 
the combined data with the value of e fixed at zero and found the fit confidence level to be 
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FIG. 17. cos( 



and (/) distributions of 7r"'"7r~ system in the e" frame in the exclusive 



T(25') T(lS')7r"*"7r measurement (corrected for acceptance). 



40.2%. Using the values of the from the two combined fits, we performed the F— test^for 
the significance of the (i— wave contribution. We calculate = A^VXn = 3.5/0.929 = 3.77 
for n = 21 d.f. which means that adding the (i— wave to the fitting function significantly 
improves the fit and corresponds to a 7% probability that the parent distribution does not 
have the ci— wave term. 

The spatial orientation of the vrvr system in the T(25') frame is consistent with isotropy 
(Figure which implies that there is no significant contribution from a "relative" D— wave 
(L = 2). 



VIII. TRANSITION T(25) T{lS)ri 

In our analysis of this transition we used the decay modes where the T(IS') decays into 
a lepton pair (e or fi) and the 1] decays via one of the modes: t] — >■ Svr^ — >■ 67, rj — >■ 27, 
1] — s> 7r+7r^7r° 7r+7r^77, or — 7r^7r~7 (the total branching ratio of these four modes is 
98.2%). The selection criteria common to the four 77 decay modes are: (1) requirements on 
the leptonic pair which are the same as in our exclusive T(2S') T(lS')7r7r analyses, (2) a 
requirement on the t] candidate momentum < 0.2 GeV/c, and (3) a requirement on the 
dilepton invariant mass 9.21 GeV/c^ < mu < 9.71 GeV/c^. 

For the modes 77 37r° 67 and — » 27 the following additional criteria are applied: 
(1) photon requirements as in the exclusive T(2S') T(lS')7r°7r° analysis except that the 
energy of 7's from r] — > Svr^ must satisfy < 0.33 GeV and those from 77 — ^ 27 must satisfy 

< 0.6 GeV, (2) there should be two good charged tracks in the event, (3) the number 
of showers in the calorimeter unmatched to charged tracks should be less then seven (for 



^See, for example, Bevington, P.R. "Data reduction and error analysis for the physical sciences." 
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FIG. 18. Signal from T(2S) T{lS)rj,T{lS) — e+e in different ij decay modes. 
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FIG. 19. Signal from T(25) ^ T(15)r?,T(15) ^ in different r? decay modes. 



22 



rj Svr") or three (for rj 27), (4) for rj —>■ Stt'^ the 7r° candidate momentum must satisfy 
Pt,o < 0.3 GeV/c, and (5) for 77 — 27 the cosine of the angle between the two photons must 
satisfy cos9^^ < —0.85 to reduce the background from the QED process e~^e~ — >• 77e"'"e~ 
(since the r^'s are produced almost at rest, the daughter 7's are close to being back to back). 

In the modes t] tt+tt^tt^ 7r+7r^77 and t] — > 7r+7r^7 we require: (1) the charged pions 
must pass the same cuts as in the exclusive T{2S) — > T(lS')7r+7r^ measurement, (2) there 
must be four good charged tracks in the event, (3) there must be three (for r] ti^ti^tt^ 
7r+7r~77) or two (for rj — »• 7r"'"7r~7) showers in the calorimeter unmatched to charged tracks, 
and (4) the cosine of the opening angle between charged pions must satisfy cos9.„+-„- < 0.9 
to suppress background from QED processes with gamma conversion 7 —>■ e^e~ where the 
e"'"e~-pair fakes a Tr+vr^-pair. 

We look for a signal in t he scatter plots of t he invariant mass of the i] candidate vs. 



the missing mass Mmiss = y{Mr{2S) ~ Ey^Y ~ p'^ which are presented in Figure |18| for the 
ee channel and in Figure |19| for the /i/i channel (the boxes denote our signal regions which 
are optimized using a Monte Carlo simulation). In Table ^ we list the number of observed 
events for the decay channels under consideration along with the detection efficiencies of 
each individual channel determined from Monte Carlo simulation. 



TABLE XI. Numbers of observed events and efficiencies for the T{2S) T{lS)r],T{lS) — > 
measurement in different r] decay modes. 



ee channel n/j, channel 



Decay mode 


BR 


jyobserved 


Efficiency (%) 


j\^observed 


Efficiency (%) 


r/ 37rO 


0.319 





2.4 ±0.3 





2.3 ±0.3 


7/ ^ 27 


0.389 


13 


38.4 ± 1.6 


3 


46.9 ± 1.9 


r] vr+TT^vr'^ 


0.236 





8.9 ±0.8 





10.5 ±0.9 


rj — > 7r"'"7r~7 


0.049 


1 


17.5 ±2.0 


2 


22.0 ±2.2 



To convert the numbers from Table ^ into branching ratios or upper limits we have 
to consider the sources of possible background contamination. The T(2S') — > T(lS')7r+7r~ 
transition with initial or final state radiation can mask the T(2S') — > T(1S')?7 transition with 
7] — *• 7r"'"7r~7 and the transition T(2S') T(15')7r'^7r'^ where two photons from different vr^'s 
escape detection can mask the r] transition with r] — > 77. To estimate these two backgrounds, 
we used our Monte Carlo sample of the exclusive dipion transitions with 50,000 events 
in the charged mode and 40,000 events in the neutral mode which we subject to the t] 
transition selection criteria. After scaling, we found the background to be 0.2(0.2) events in 
the rj — > 7r+7r~7 ee(/i/i) channel and 0.3(0.6) events in the 77 — 77 ee(/i/i) channel. We did 
not observe any background events in the rj —>■ Sn^ or rj ^ tt+tt^tt'' channels. Another source 
of possible background are the cascade radiative decays T(2S') 7x5 -y-yT^lS). This 
contamination was estimated based on a 15,000 Monte Carlo sample of the cascade radiative 
decays. We found no background events from this source. To estimate the background from 
radiative QED and other possible nonresonant processes we used a data sample corresponding 
to an integrated luminosity of 5.17 pb~^ of e~^e~ annihilations taken at = 9.98 GeV, just 
below the T(2S') resonance. After scaling for luminosity and energy differences we found 
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14.2(0.0) background events for the rj —>■ 'yj mode in the ee{fifi) channel and no background 
events for the three remaining 77 decay modes. The resuhs of the background study are 
summarized in Table pCll[ 

TABLE XII. Number of expected background events for the T(25) T{lS)ri,T{lS) l+l' 
transition in different ij decay modes for our T(2S') resonance data sample. 



Decay mode 




Sources of background, events in ee(/x/i) channel 


Total 






77 cascade 


T(25), continuum 


r/ SttO 


0(0) 


0(0) 


0(0) 


0(0) 


0(0) 


r/ ^ 27 


0(0) 


0.2(0.2) 


0(0) 


14.2(0) 


14.5(0.2) 


r] — > 7r+7r~7r'' 


0(0) 


0(0) 


0(0) 


0(0) 


0(0) 


r] — >■ 7r"'"7r~7 


0.3(0.6) 


0(0) 


0(0) 


0(0) 


0.3(0.6) 



Although the above study shows that in the /i/i channel the expected number of events 
from background processes in the signal region is not consistent with the number of observed 
events, some of the signal events lie very close to the signal box boundary which leads us 
to interpret our signal candidates as smearing of background events into the signal region. 
Therefore we (conservatively) do not calculate a branching ratio but set an upper limit. 
Because the mode 77 ^ 77 in the ee channel is so "noisy" we exclude it from the procedure 
of setting the upper limit on the branching ratio. After taking into account the errors on 
efficiencies, on the T(IS') leptonic and t] branching ratios, we set the following upper limit:Q 

B{T{2S) T{1S)7]) < 0.0028 (90%aL.) 

The results from other experiments are given in Table PCII4 

TABLE XIII. Upper limits on B{T{2S) T(15)r?) (90%aL.). 



CLEO g < 0.010 

Crystal Ball |2| < 0.007 

ARGUS ig < 0.005 

CUSB ig < 0.002 

this analysis < 0.0028 

In the multipole expansion of the gluon color field, tttt transitions proceed via El ■ El 
emission but the lowest order transition allowed by the quantum numbers of the ?7-meson 
is El ■ M2 or Ml ■ Ml emission. This results in a suppression of the rj transition com- 



pared to the vr^vr" transition by a factor of ~ 5 ■ 10~ ||2^, so the branching ratio for 
T(2S') T{lS)ri is expected to be around 0.001, below the current upper limit. Since 
for the chromomagnetic transitions the transition amplitude varies as rn~^^^i^, the ratio 



^''To calculate an upper limit on the number of signal events we follow the procedure suggested 
by PDG and include the systematic errors according to [p^]. 
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B{T{2S) T{lS)ri)/B{T{2S) T(l^)7r+7r-) should be substantially smaller than the 
ratio B{ip{3685) V^r/)/i3(i'(3685) ^tt+tt-) = 0.083. Yan @] obtained the formula: 

^ r(T(2g) ^ T{lS)r]) ^ /m^Y /prV ^ J_ 
^^'^ r(V'(3685) ^ ^/-T^) ~ VmJ ~ 275 

where px and are the decay momenta. Our experimental value is r^jc < 1/61? using 
rtot(V'(3685)) = 277 keV and B{i/j{3685) i/jt]) = 0.027; this is 15 times smaller than 
the suppression expected from phase space alone (a factor of four). Our results are clearly 
consistent with the multipole expansion formalism. 



IX. TRANSITION T(2S') T(15)7r0 

We also studied the isospin violating transition T(25') — > T(lS')7r° with T(15') — > 
and 7r° — *• 77. The same set of selection criteria as in the exclusive T{2S) —>■ T{lS)7i7i study 
is applied to lepton candidates and the same set of selection criteria on photons that was 
used for direct reconstruction of r/'s from two 7's in the T{2S) T{lS)ri study is applied 
here. Additional requirements are: (1) p^o < 0.6 GeV/c, (2) the number of good charged 
tracks equals two, (3) the number of showers unmatched to tracks less than three, (4) the 
cosine of the angle between the vr" and the dilepton system cos9t,ii < —0.9 (to reduce the 
background from QED processes), and (5) 9.21 GeV/c^ < mu < 9.71 GeV/c^ where mu is 
the dilepton invariant mass. 

As in the search for the t] transition, we look for a signal in the scatter plot of the vr^ 



invariant mass m^o vs. the missing mass Mmiss = y (^t(25) — E^^oY — p^o- In Figure |2^ the 
scatter plots from the T{2S) resonance data sample are displayed for both dilepton channels 
(Monte Carlo simulation is used to optimize the signal regions denoted by the solid boxes). 



TABLE XIV. Number 


of observed events 


and efficiencies for the T(25) T(15')7r^, 


T(IS') transition. 








jyobserved 


Efficiency (%) 


ee 


9 


29.3 ±0.8 




6 


36.3 ±0.9 



Within the signal region, we find 9 events in the ee channel and 6 events in the /i/x 
channel. The efficiencies, which are based on Monte Carlo simulation, are given in Table ^TV. 

We used a "grand side-band" technique, to estimate the background: we counted the 
events in the "grand side-band" (in Figure |20| it is all the area outside the signal box for 
the ee channel and a vertical strip between 9.41 GeV and 9.51 GeV in Mmiss, excluding the 
signal box, for the fifi channel) and extrapolated the background event yield into the signal 
region. The results are given in Table pCV. 



As seen in the table, using the "grand side-band" subtraction technique we expect 12.9 
background events compared to the total of 15 observed events. This corresponds to an 
upper limit: 
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FIG. 20. Scatter plot of vr'' invariant mass vs. missing mass for T{2S) T(lS')7r'', 
T{1S) from T{2S) resonance data. 

TABLE XV. Numbers of the events from the "grand side-band" subtraction technique. 

background 
siQnal—reqion 



N. 



ee 



sideband 

37 



N. 



signal— region 

"9 
6 



Estimated A^, 



8.4 
4.5 



B{T{2S) ^ T(15)7r0) < 0.0011 (90%C.L.) 
This is the lowest upper hmit on the 7r° transition to date. The only other experiment that 



studied this transition was Crystal Ball (Table |XV|). 

TABLE XVL Upper Hmits on B{T{2S) T{IS)tt^) {MC.L.). 



Crystal Ball l2| 


< 0.008 


this analysis 


< 0.0011 



The T{2S) —>■ T(15')7r'^ transition can occur because of a breaking of the isotopic sym- 
metry due to the mass difference between the u and d quarks, and its rate is expected to be 
lower than the T[2S) —>■ T{lS)ri rate. In the context of the multipole expansion, this ratio 
is given by p8| : 

r((2g) ^ (lg)7r°) ^ / m,-m„ y fmA' (pA' 
'■'^"/''^ Tii2S) ^ {lS)v) [m, + mj [mj [pj 

With {rriii — mu)/{ma + m„) ^ 0.3 this gives r^o/r, ~ 0.022 for charmonium which is 
in reasonable agreement with the experimental value of 0.037. For bottomonium we have 
r^o/^ ^ 0.14 and F(T(2^) ^ T(l^)7rO) ^ 0.003 keV (using F(T(2^) ^ T(l^)r/) = 0.02 keV 
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from Kuang-Yan [^) which is more than an order of magnitude below our upper hmit of 
0.048 keV. 



X. SUMMARY 



We have measured various experimental quantities for the hadronic transitions from the 
T{2S) to T(IS') including branching ratios, the dipion invariant mass spectra, and angular 



distributions. Using the PDG value for the full width of the T{2S) resonance F = 44 keV [30 



we also calculate the partial widths for the corresponding transitions. Table PCVII] reports 



our measurements of the branching ratios and partial widths compared with previous world 
averages and theoretical calculations by Kuang and Yan . Our results are consistent with 
previous experiments as well as theoretical predictions. We determine an upper limit on the 
branching ratio of T{2S) T{lS)ri and set a new upper limit on the branching ratio of the 
T{2S) T{1S)tt^ transition. 

TABLE XVII. Summary of the branching ratios and rates of hadronic transitions of T(25). 







Branching 


ratio 


Rates (keV) 


Decay 




Experiment 


World Avg. 


Experiment 


Kuang-Yan 


T{2S) - 


T{lS)7r+TT- 


0.192 ± 0.002 ± 0.010'' 


0.185 ±0.008 


8.4 ±0.5 


8.8 


T{2S) - 




0.092 ± 0.006 ± 0.008 


0.088 ±0.011 


4.0 ±0.4 


4.4 


T{2S) - 




< 0.0028 


< 0.002 


< 0.12 


0.02 


T{2S) - 


T(15)7r0 


< 0.0011 


< 0.008 


< 0.048 


0.003 



" average over the exclusive and inclusive measurements. 



We also calculate the leptonic branching ratios of the T(IS'): B^e = B(T{1S) e+e ) = 
0.0229 ± 0.0008 ± 0.0011 and B^^ = B{T{1S) fi+fx') = 0.0249 ± 0.0008 ± 0.0013 which 
are in good agreement with PDG values. 

The dipion invariant mass spectrum we observe in T(25') T(lS')7r7r transitions is well 
described by both the model of multipole expansion of the gluon color field and the model 
of Novikov, Shifman, Voloshin and Zakharov who used the general form of the QCD field 
tensor G^^, to obtain a hadronization matrix element in the chiral limit P,p^. 

The angular distributions of the final state particles in T(2S') — > T(lS')7r7r show a strong 
s— wave dominance, as expected from theory. A small ci— wave contribution on the order of 
4% may be present in our data. 
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